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ABSTRACT 

We present high spatial resolution ground-based broad-band imaging, Ha+[NII] 
narrow-band imaging and long-slit spectroscopy for the E4 radio galaxy NGC 7052, 
which has a nuclear dust disc. We detect ionized gas with a LINER spectrum, residing 
also in a nuclear disc. The gas rotates rapidly (■Urot ~ 250 km s"-*^ at 1.5"). The emis- 
sion line widths increase towards the nucleus (nuclear FWHM 500 — 600 km s"-*^). 

The images are well fit by an axisymmetric model with the gas and dust in a disc of 
~1.5" radius (~340 pc) in the equatorial plane of the stellar body, viewed at inclination 
i 70°. The flux distribution of the ionized gas is very centrally concentrated, and 
is unresolved inside the seeing radius. The central cusp-steepness a of the stellar 
luminosity density [j oc for t" — > 0) is not well constrained by the observed surface 
brightness distribution, because of the dust absorption. 

We model the gas kinematics with the gas on circular orbits in the equatorial plane, 
with a local velocity dispersion due to turbulence (or otherwise non-gravitational mo- 
tion). The circular velocity is calculated from the combined gravitational potential of 
the stars and a possible nuclear black hole. The observed gas rotation curve is well fit 
by a model with a = —1.3, either with or without a black hole. Turbulent velocities 
> 300 km s"-*^ must be present at radii < 0.5" to fit the observed nuclear line widths. 
Seeing convolution of a Keplerian rotation curve around a 10^ Mq black hole can fit 
the observed widths without turbulence, but such a model predicts nuclear emission 
line shapes with pronounced peaks at v ±300 km s"-*^, which are not observed. Mod- 
els with both a black hole and gas turbulence can fit the data well, but the black hole 
is not required by the data, and if present, its mass must be < 5 x 10® M©. 

Although this upper limit is not very stringent, it is already ~5 times smaller than 
the black hole mass inferred for M87 from HST data. We show that HST observations 
of NGC 7052 should improve significantly the constraints on the mass of any possible 
black hole. Kinematic observations of nuclear gas discs are likely to become a widely 
used tool in the search for massive black holes in galactic nuclei. The modelling and 
analysis techniques presented here will be useful for the interpretation of such data. 

Key words: galaxies: active - galaxies: elliptical - galaxies: individual: NGC 7052 
- galaxies: kinematics and dynamics - galaxies: nuclei - line: profiles. 



1 INTRODUCTION 



clear-cut for none of them, due to finite spatial resolution 
and lack of knowledge on the shape of the stellar orbits. 
An alternative is to study the dynamics of ionized gas, if 
present. This can be very complicated if the gas morphology 
is irregular or unresolved, or if there are significant contri- 
butions from, e.g., inflow or outflow. This is often the case 
for the narrow-line and broad-line regions of active galac- 
tic nuclei. However, if there is observational evidence that 
the gas moves on (approximately) circular orbits in a disc, 
the analysis is straightforward and less complicated than for 
the stellar kinematics. M87 provides a good example. Van 



It is generally believed that active galaxies are powered by 
the presence of a massive central black hole (e.g., Rees 1984). 
In fact, black holes are expected to be present in the centres 
of many quiescent galaxies as well (e.g., Chokshi & Turner 
1992). However, direct dynamical evidence for the pres- 
ence of black holes in individual galaxies is scarce. Stel- 
lar kinematic studies (e.g., Kormendy 1993; van der Marel 
et al. 1994a, b,c) have provided tentative evidence for black 
holes in a handful of nearby galaxies. However, the case is 
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der Marel (1994c) detected high velocities near the nucleus 
which he modelled as gravitational motion around a black 
hole. Subsequent observations with the Hubble Space Tele- 
scope (HST) by Ford et al. (1994) and Harms et al. (1994) 
resolved the gas spatially, and showed it to be in a rapidly 
rotating disc. The observed gas velocities in M87 provide 
the currently best evidence for a black hole in any galactic 
nucleus (with the possible exception of our own galaxy). 

NGC 7052 is an E4 galaxy with apparent magnitude 
ms = 13.4 (de Vaucouleurs et al. 1991), at heliocentric ve- 
locity ■Uhei = 4687 km s"'' (Faber et al. 1989). It is a known 
radio source with a core and jet, but no lobes (Morganti 
et al. 1987). The radio axis is offset from the photometric 
major axis of the galaxy by ~50° (MoUenhoff, Hummel & 
Bender 1992). It is also an infrared source. It was detected 
in all four IRAS bands and has a far-infrared luminosity of 
~10^° Lq (Impey, Wynn- Williams & Becklin 1990). Opti- 
cal images show absorption by dust near the nucleus. Very 
high spatial resolution (seeing FWHM 0.33") CFHT images 
of the nucleus were presented by Nieto et al. (1990). These 
show that the dust resides in a lane aligned with the ma- 
jor axis (see also Figure 1), suggesting it to be in a nearly 
edge-on disc with ~1.6" radius. Molecular gas was detected 
in NGC 7052 using the CO (1-0) line (Wang et al. 1992). 
Atomic gas has not been detected (Mhi < 5 x 10* Mq; Dres- 
sel, Bania & O'Connell 1982). Morganti, Ulrich & Tad- 
hunter (1992) observed NGC 7052 through a narrow-band 
filter centred on Ha and the adjacent [Nil] lines. Ionized 
gas was detected around the nucleus, but its morphology 
was unresolved at the ~2" resolution of their observations. 

The presence of a nuclear dust disc in NGC 7052 sug- 
gests that the ionized gas might also reside in a disc, since 
the occurrence of dust and gas in elliptical galaxies are often 
closely coupled (Bertola 1992; Buson et al. 1993; Goudfrooij 
et al. 1994). If so, observations of the gas kinematics could 
yield useful constraints on the nuclear mass distribution. 
For this reason we observed NGC 7052 in August 1992 with 
the 4.2m William Herschel Telescope (WHT) on La Palma. 
Broad-band images (Section 2) and long-slit emission- and 
absorption-line spectra (Section 4) were obtained in excel- 
lent seeing conditions (FWHM w 0.6"). A narrow-band 
image was later also obtained (Section 3). Our observa- 
tions show that the ionized gas indeed resides in a rotating 
disc. Interpretation of the observed motions requires care- 
ful modelling to correct for the effects of seeing and dust 
absorption. In Section 5 we model the stellar light distri- 
bution, the geometry of the ionized gas, and the geometry 
and optical depth of the dust. We thus fit the broad- and 
narrow-band images and the integrated emission line fluxes 
derived from the spectra. In Section 6 we model the ob- 
served gas kinematics. We show that a central black hole 
is not required to fit the data, but that some form of gas 
turbulence is. We put an upper limit on the mass of any 
possible black hole (or other dark compact nuclear mass). 
For our best-fitting models we calculate predictions for fu- 
ture spectroscopic observations with the HST. Conclusions 
are presented in Section 7. 

Values quoted throughout this paper are based on a 
Hubble constant Ho = 100 /i km s"^ Mpc"^ with h = 1.0. 
The particular choice adopted for h does not influence the 
predictions of any of our models. It merely sets the length, 
mass and luminosity scale of the models in physical units. 



Figure 1. ii-band image of tlie central region of NGC 7052, 
obtained with seeing FWHM S = 0.67". The directions of North 
and East are as indicated. There is a dust lane in the inner 2", 
aligned with the major axis of the galaxy. 

Specifically, distance, length scales and masses scale as , 
while mass-to-light ratios scale as h. 

2 BROAD-BAND IMAGING 

2.1 Data acquisition and reduction 

Three il-band images of NGC 7052 were obtained with the 
GEC5 CCD at the Cassegrain focus Auxiliary Port of the 
WHT. The CCD pixels are O.IO" x O.IO", with a 38" x 58" 
field of view. The image reduction was done using standard 
IRAF tasks. 

The seeing point spread functions (PSFs) of the obser- 
vations were determined from bright stars on the images. 
The PSF shape could be well described by a sum of two 
Gaussians: 

2 /r, 2 2 /r, 2 

PSF(7-) = ^le"" -^^'i +^26-'' (1) 

with fixed ratios of the dispersions, 1T2/1J1 = 1.92, and am- 
plitudes A2/A1 = 0.29, of the two components. The PSF is 
normalized for Ai = 0.0774:3/ af , and is fully specified by its 
FWHM S = 2.644(71. 

Our best image is a 600 sec. exposure with seeing 
FWHM S = 0.67" ± 0.02". The analysis below is based 
on this image. The other two images were not used be- 
cause of poor signal-to-noise ratio [S/N) and poor seeing, 
respectively. Figure 1 shows the central region of NGC 7052. 
There is a sharp dust lane aligned with the major axis of the 
galaxy, in agreement with the (higher spatial resolution) im- 
age of Nieto et al. (1990). 

The WHT autoguider, which was locked onto a bright 
star during each exposure, provides an independent mea- 
sure of the seeing FWHM, 5auto. It does not properly take 
the full PSF shape (1) into account though, and its esti- 
mates must be calibrated. Our images imply 5/5auto ~ 0.76, 
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Figure 2. i?-band surface brightness, ellipticity, major axis position angle and cos 46 Fourier coefficients of the isophote deviations from 
ellipses, determined from the image in Figure 1. The cos 36, sin 36 and sin 46 coefficients were found to be close to zero, and are not 
shown. The abscissa is the radius along the major axis, measured with respect to the centre of the outer isophotes. The observed surface 
brightness at radii < 2.5" cannot be meaningfully approximated by ellipses. The stellar surface brightness and dust geometry in this 
nuclear region are modelled in Section 5. 



in good agreement with the calibration derived by van der 
Marel et al. (1994a). Calibrated autoguider estimates are 
used in Section 4 to estimate the seeing FWHM for our 
WHT spectra. 



2.2 Isophotal Analysis 

We analysed the isophote shapes with the ellipse fitting soft- 
ware of Jedrzejewski (1987), in the manner described by van 
den Bosch et al. (1994). Figure 2 shows the il-band surface 
brightness, ellipticity (e), major axis position angle (PA) and 
cos 45 Fourier coefficients of the isophote deviations from el- 
lipses, as functions of radius along the major axis. 

Accurate isophotal parameters could not be determined 
at radii > lO", because the galaxy was offset towards one 
corner of the CCD. This allows an accurate sky subtraction 
(using the measured brightness in the opposite corner of the 
CCD) at the cost of a more limited radial extent of the data. 
No isophotal results are shown for radii < 2.5", where the 
isophotes cannot be meaningfully approximated by ellipses. 
Models for the stellar surface brightness at small radii, and 
for the geometry and optical depth of the dust are discussed 
in Section 5. 

Between r = 2.5" and r = 10" the ellipticity of the 
isophotes increases from ~0.3 to ~0.4. The position an- 
gle is constant and the higher-order Fourier terms are close 
to zero. These results are consistent with previous studies 
by Bender, Dobereiner & MoUenhoff (1988) and Nieto et 
al. (1991). Both find that the ellipticity increases further to 
e Ri 0.55 at ?• = 50". There is no agreement in the liter- 
ature on the isophote deviations from ellipses for r > 20". 
Bender et al. find the isophotes to be boxy (negative cos 45 
coefficient), Nieto et al. find the isophote deviations to be 
close to zero, while Ebneter, Djorgovski & Davis (1988) ar- 
gue that NGC 7052 might be an SO (which would seem to 
be consistent with the increasing ellipticity with radius; see 
the arguments in van den Bosch et al. 1994). 



3 NARROW-BAND IMAGING 

Marcella CaroUo kindly obtained narrow-band images of 
NGC 7052 for us in July 1994, using the 1.5m Danish Tele- 
scope at ESQ, La Silla. The CCD pixels are 0.377" x 0.377", 
with a 6.4' x 6.4' field of view. Two 30 min. exposures were 
taken, one 'on-band' and one 'off-band'. The on-band filter 
covered the region around Ha (Arest = 6562. 68A) and the 
adjacent [Nil] lines (Arest = 6548.03, 6583. 4lA), at the red- 
shift (AA w IO2A) of NGC 7052. The off-band filter was 
centred on A w 6550A. Both filters have an effective width 
of ~60A. 

After bias, dark current, cosmic ray and sky subtrac- 
tion, and flat-fielding, the frames were shifted to a common 
coordinate system using stars as reference points. No ro- 
tation of the frames was necessary. After measuring the 
seeing FWHM from stellar images, the best seeing frame 
was convolved with a Gaussian to yield an 'effective' see- 
ing FWHM equal to the seeing FWHM for the worst seeing 
frame (5 = 1.41"). Subsequently, the images were scaled to 
the same surface brightness at radii > lO" (were there is no 
line emission according to our spectra in Section 4.1), and 
the off-band image was subtracted from the on-band image. 

The resulting Ha+[NII] image is shown in Figure 3. 
The innermost gas contours are flatter than the stellar con- 
tours, in spite of the relatively poor seeing of the narrow- 
band image (which tends to make the gas contours rounder). 
Hence, the ionized gas most likely resides in a disc. The ma- 
jor axes of the stellar and gas components are aligned within 
2°. The centre of the Ha + [Nil] emission is offset from the 
centre of the outer isophotes of the galaxy by ~ 0.4". It 
is shown below (Section 5) that this can well explained as 
due to absorption of part of the line emission by the dust, 
which itself is not entirely symmetrically distributed with 
respect to the centre. There might be an additional con- 
tribution from the fact that the ionized gas rotates rapidly 
(Section 4.2), causing the transmission of the emission line 
flux by the on-band filter to vary with position. 
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spectra along the slit with —5" < < 5" were summed. 
This spectrum was then analysed with the pixel space fit- 
ting method of van der Marel (1994c), using the spectrum of 
a K5 III star obtained with the same instrumental setup as 
template. This yields an 'average' stellar velocity dispersion 
for the central region of NGC 7052: = 266 ± 26kms"\ 
This is consistent with the higher quality data of Wagner, 
Bender & MoUenhoff (1988), who inferred a central velocity 
dispersion of ~270 km s"'' . They found that the stellar ro- 
tation velocities are comparable along the major and minor 
axes, indicative of triaxiality, but that the overall contribu- 
tion of rotation to the dynamical support of the system is 
negligible (|Kot/o-| < 0.2). Unfortunately, the 5/JV of our 
absorption-line spectrum is too low to study the variation 
of the stellar kinematics as a function of radius. 



Figure 3. Map of the Hct+[NII] emission (thick contours) over- 
plotted on the observed broad-band surface brightness distribu- 
tion from Figure 1 (thin contours). The narrow-band image was 
obtained in relatively poor seeing (FWHM S = 1.41"). The in- 
nermost gas contours are none the less flatter than the contours 
of the stellar component, indicating that the gas most likely re- 
sides in a disc. The offset between the centre of the HcK-t-[Nll] 
emission and the centre of the outer stellar isophotes can be well 
understood as due to absorption of part of the line emission by a 
slightly asymmetric dust distribution (cf. Figures 8 and 10 below). 



4 LONG-SLIT SPECTRA 

Long-slit spectra of NGC 7052 were obtained with the 
R1200R grating and the EEV3 CCD on the red arm of 
the ISIS spectrograph of the WHT. The dispersion is 0.38 
A/pixel, the spatial scale 0.33"/pixel. A slit width of 
0.8" was used, yielding an instrumental velocity dispersion 
finstr ~ lOkms"''. Exposures of arc lamps were taken be- 
fore and after each galaxy exposure to allow accurate wave- 
length calibration. The galaxy was centred on the slit at 
the beginning of each exposure using a TV camera with an 
il-band filter. The autoguider was equipped with an /-band 
filter. Differential atmospheric refraction between the R and 
/-band is negligible (< O.l"). The data were reduced as in 
van der Marel et al. (1994a). Standard steps included: bias 
subtraction, flat-fielding, cosmic ray removal, wavelength 
calibration, logarithmic wavelength rebinning, rebinning to 
correct for S-distortion, sky subtraction and dark current 
subtraction. The wavelength dependence of the efficiency 
of the telescope-detector combination was corrected using 
observations of a spectrophotometric standard star. 

4.1 Absorption-line spectra: stellar kinematics 

One 50 min. major axis (PA = 65°) absorption-line spec- 
trum was obtained from 8390A to 8820A, covering the Ca 
II IR triplet region. The seeing FWHM was 5 = 0.84" 
(obtained from the autoguider estimate; see Section 2). To 
obtain sufficient S/N for a stellar kinematic analysis, all 



4.2 Emission-line spectra: ionized gas kinematics 

One 30 min. major axis (PA = 65°) and one 30 min. minor 
axis (PA = 155°) spectrum were obtained from 6435A to 
6880A, around IIa+[NII]. The seeing FWHM was S = 0.57" 
for the major axis spectrum and S = 0.84" for the minor axis 
spectrum (obtained from the autoguider estimates). Hence, 
the major axis spectrum was taken under very good seeing 
conditions, although not exceptional for La Palma (see, e.g., 
Zeilinger, M0ller & Stiavelli 1993). 

We performed the following test to verify the seeing es- 
timate for our major axis spectrum. The /l-band image in 
Figure 1 was 'overlaid' with a hypothetical slit along the ma- 
jor axis, with the same width and pixel size as for our spec- 
troscopic observations. Upon summing the /l-band intensity 
over each pixel, this yields the predicted intensity profile 
along the slit for spectroscopic observations with S = 0.67" 
(the seeing FWHM for the /l-band image). This predicted 
profile was found to be less centrally peaked than the ob- 
served intensity profile along the slit, consistent with the 
fact that the major axis spectrum had S < 0.67". 

To study the emission lines one needs to subtract the 
spectrum of the underlying stellar component. Fitting this 
component merely with a low-order polynomial can intro- 
duce systematic errors (e.g., Goudfrooij et al. 1994). Hence, 
we model the stellar component as the convolution of a tem- 
plate spectrum (the spectrum of a K3 III star obtained with 
the same instrumental setup) with a Gaussian velocity dis- 
tribution. The latter was taken to be the same as determined 
in Section 4.1. This is adequate, because the stellar velocity 
dispersion does not vary much, and the stellar rotation ve- 
locities are small, over the spatial region where the emission 
lines are observed (Wagner et al. 1988). 

No emission lines were detected outside the central 2", 
in agreement with Figure 3. Figure 4 shows the emission- 
line spectra for three representative radii along the major 
axis. There are no residual absorption lines, indicating 
that the starlight has been properly subtracted. In addi- 
tion to Ha and the [Nil] lines, we detect the [SII] lines 
with Arest = 6717.47, 6730. 85A. The nuclear line ratios, 
[NII]6583/Ha w 1.7 and [SII]67i7+673i /Ha w 1.0 are char- 
acteristic for a low ionization emission-line region (LINER), 
and are significantly larger than those seen in HII region-like 
objects (Veilleux & Osterbrock 1987). 

The peak [Nil] 6583A line flux in our nuclear spectrum 
is 1.8 X 10"^^ ergcm"^ s"^ A"\ The integrated Ha+[NII] 
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Figure 4. Emission line spectra at three representative radii 
along the major axis (0.61", —0.05", —0.71"), after subtraction 
of the stellar contribution to the spectrum. The spectra are char- 
acteristic for a LINER. The vertical scale is in arbitrary units. 

flux summed over all radii along the major axis is 1.8 x 
10"''^ erg cm~^ . Morganti et al. (1992) measured a total 
Ha+[NII] flux for NGC 7052 of 2.4 x 10"^^ erg cm"^ s"^ . 

To determine the ionized gas kinematics we study the 
Ha and [Nil] emission lines. The [SII] lines are not used, 
because of their lower S/N. We obtain smooth fits to the 
Ha and [Nil] lines using a two step process. In each step a 
suitably parametrized function is fit to the data using the 
Levenberg-Marquardt method (e.g.. Press et al. 1992). The 
quantity that is minimized takes read-out and photon 
noise into account. In the first step all three lines are as- 
sumed to be Gaussian, and the three Gaussians are fit to 
the data simultaneously. The ratio of the amplitudes of 
the [Nil] lines is fixed to 3, since this is the ratio of their 
transition probabilities (Osterbrock 1989). Since both [Nil] 
lines originate from the same gas component, their velocity 
dispersions are constrained to be equal, and their velocity 
difference is fixed to Av = cAlnA = 1615.5 km s"'^ . In the 
second step each line is fit separately using a function that is 
more general than a single Gaussian, namely the sum of two 
Gaussians. This has no physical significance, but is merely 
a convenient fitting function. Before fitting either line the 
other two lines are subtracted, using the Gaussian fits ob- 
tained in the first step. This procedure works well since 
there is only little overlap between different lines (Figure 4). 
Figure 5 shows the best double Gaussian fits to the major 
axis Ha lines after subtraction of Gaussian fits to the [Nil] 
lines, to illustrate the procedure. 

From the best-fitting double Gaussian to each line we 
determined the velocity moments up to order 2 (the in- 
tegrated line flux, mean velocity and velocity dispersion) 
and their formal errors. These are convenient quantities for 
comparison to the predictions of theoretical models. (Al- 
ternatively, the width of the emission lines could have been 
characterized by their FWHM, rather than by their disper- 
sion IT. For a Gaussian profile these are related by FWHM 
= VS In 2(7.) Figure 6 sh ows the results for the major and 
minor axes. The two individual [Nil] lines yielded consis- 
tent results. The results shown in Figure 6 are for the 
line at Arest = 6583. 4lA, which has the highest S/N. The 
ionized gas rotates rapidly along the major axis, reaching 
~250kms~'' at ?• = 1.5". The velocity dispersion of the 



Figure 5. Observed Hck emission lines at different radii along 
the major axis (thin curves), after subtraction of Gaussian fits 
to the neighbouring [Nil] lines. Thick curves are the best-fitting 
double Gaussians. The shift of the lines with radius indicates 
rapid rotation of the gas. The velocity dispersion of the gas peaks 
in the centre. 

ionized gas has a central peak. There are some differences 
between the results for Ha and [Nil]. The flux distribution 
of [Nil] is more centrally peaked than that of Ha. Along the 
major axis, the rotation velocities of [Nil] are lower and the 
velocity dispersions higher than those of Ha. 
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Figure 6. Integrated flux (in arbitrary units), rotation velocity and velocity dispersion of the ionized gas along the major and minor 
axes of NGC 7052. Results for Hck and [Nil] are indicated with different symbols. Error bars are shown for Hck only, to avoid confusion. 
The errors on the results for [Nil] are of similar magnitude. 



5 MODELLING THE STELLAR LIGHT 
DISTRIBUTION AND GAS AND DUST 
GEOMETRY 

5.1 Model assumptions 

We assume for the three-dimensional stellar luminosity den- 
sity j: (i) that it is oblate axisymmetric; (ii) that the iso- 
luminosity spheroids have constant flattening Q as function 
of radius; and (iii) that j can be parametrized as: 

j{R,z) = jo{m/b)^[l + {m/bff, m' = R'+z'Q-\ (2) 

where {R, z) are the usual cylindrical coordinates, and jo, 
b, a and /3 are free parameters. When viewed at inclina- 
tion angle i in the absence of dust absorption, the projected 
intensity contours are ellipses with axial ratio Q' , with 

Q'^ = cos^ i -\- sin^ i. (3) 

Models with the same Q' but different i have the same pro- 
jected intensity apart from a normalization constant. The 
projected intensity must generally be calculated numerically. 
Simple analytical expressions exist only for special combina- 
tions of the model parameters. 

The assumption of axisymmetry for the luminosity den- 
sity cannot be entirely correct, given that stellar streaming 



is observed along both the major and minor axes (Wagner et 
al. 1988). However, there are several indications that NGC 
7052 might in fact be close to axisymmetric, or at least only 
mildly triaxial. There is no isophotal twisting (Figure 2), 
the dust (Figure 1) and ionized gas (Figure 3) are aligned 
with the major axis of the galaxy, and there is no significant 
rotation of the gas along the minor axis (Figure 6). 

The assumption of constant flattening for the luminos- 
ity density was adopted because it simplifies the dynam- 
ical modelling. However, it does not fit the observations 
very well. The observed ellipticity of the isophotes increases 
steadily with radius beyond ~2.5" (Figure 2). None the less, 
we consider our assumption to be adequate. First, the main 
thing we want to include in the model is that NGC 7052 
is flattened; changes in flattening are only a second order 
effect. Second, what is most important for our dynami- 
cal modelling is the flattening in the central ~2.5", which is 
poorly constrained observationally because of the absorption 
by dust. Adopting models with variable flattening would 
not resolve this degeneracy. Third, the influence of ellip- 
ticity changes can be assessed by comparing the dynamical 
predictions of models with different constant flattening. We 
did this, and found no significant change in any of the major 
conclusions of our paper. In the following we adopt Q' = 0.7 
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Figure 7. The top panels show the observed minor axis i?-band surface brightness profile in the direction least affected by dust 
absorption, as function of distance to the centre of the outer isophotes. The solid curves are fits of the seeing convolved projected 
intensity of our models to the data with r > l". The dashed curves show the predictions before seeing convolution. Two representative 
choices for the cusp steepness ck are shown, ck = (a constant surface brightness core) and ck = —1.3 (a central power law cusp). The 
lower panels show the residuals between the model and the observations, as resulting from dust absorption. The cusped model requires 
a larger optical depth for the dust, and a dust distribution extending to a larger radius. 



(i.e., £ = 0.3), which fits the observed ellipticity at ?• = 2.5". 

The parametrization (2) for the luminosity density has 
no particular physical significance. It was chosen because it 
reproduces the generic behaviour of the surface brightness 
profiles of elliptical galaxies observed with the HST: a cen- 
tral power law cusp in the inner parts with a smooth change 
at some intermediate radius to a steeper power law fall-off 
at larger radii (e.g.. Crane et al. 1993; Ferrarese et al. 1994; 
Kormendy et al. 1994). Models with a = have a constant 
luminosity core. Models with a > — 1 have a constant sur- 
face brightness core when projected onto the sky. Models 
with a < — 1 have a central power-law surface brightness 
cusp when projected onto the sky. 

For lack of information on the distributions of the ion- 
ized gas and dust along the line of sight, we assume both 
components to be infinitesimally thin, residing in the equa- 
torial plane of the galaxy. Furthermore, we assume the ion- 
ized gas flux F to be circularly symmetric, F = f (-R). We 
consider distributions of the form: 

= I ^1 exp{-R/Ri) + h exp{-R/R2), {R < R) 
0, [R > R) . 

This parametrization was chosen because with a proper 
choice for the free parameters Ji , I2, Ri, R2 and R, it 
can adequately fit our data. It has no further significance. 
The flux distributions of Ha and the [Nil] lines are differ- 
ent (Figure 6), and we thus use different sets of parameters 
to describe them. The optical depth r of the dust is not 
constrained to be circularly symmetric, because it is not a 



priori clear whether the dust resides in a disc (or ring), or in 
a single filament. We determine the dust distribution from 
the data using an unparametrized technique (Section 5.3). 

The observed emission line flux fobs need not repre- 
sent the intrinsic flux F, since part of it might have been 
absorbed by dust. The fraction of the emission line flux 
that is absorbed depends on the distributions of the gas and 
dust in the direction perpendicular to the equatorial plane, 
which we do not model. We therefore include absorption 
of line emission by dust phenomenologically, by considering 
two limiting cases: fobs = f ('gas in front of the dust') and 
fobs = e~^f ('gas behind the dust'). In the latter case, r 
is assumed to be the same as in the il-band. The difference 
in extinction between the wavelengths of the il-band and 
IIa+[NII] is negligible (Osterbrock 1989). 

Our model for the gas and dust in NGC 7052 is not 
elaborate. However, it captures the essential physics of the 
problem, and as will be shown below, it allows good fits 
to be obtained to the observed dust and gas morphology, 
including its asymmetries. 

5.2 The minor axis surface brightness profile 

A useful first constraint on the parameters of the luminos- 
ity density j is obtained by fitting to the observed minor 
axis surface brightness profile in the direction least affected 
by dust absorption ('up' in Figure 1). Figure 7 shows this 
profile as function of distance to the centre of the outer 
isophotes. This centre is not the point of maximum surface 
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Figure 8. Thick contours sliow tlie (seeing-deconvolved) dust optical depth, obtained with the algorithm in Section 5.3, for i = 70° and 
two representative values of the cusp-steepness a.. The left panel shows the contours of the optical depth for ck = (contours at 0.2, 0.5, 
0.8 and 1.1), while the right panel shows the more symmetrical optical depth distribution for ol = —1.3 (contours at 0.2, 0.5, 0.9, 1.5 and 
3.0). Thin contours show the observed broad-band surface brightness distribution. 



brightness, due to dust absorption. The cusp-steepness a is 
thus not well constrained. 

The profile outside ~l" shows little sign of being influ- 
enced by dust absorption. For any fixed choice for a, the 
parameters jo, b and /3 can thus be determined by fitting 
the seeing convolved projected intensity of the model to the 
data outside Figure 7 shows the best fits for two rep- 

resentative cases: a = and a = —1.3. The residuals of 
the fit in the central arcsec are largest for the model with 
a = —1.3. This model has more light near the centre, and 
thus requires a larger optical depth for the dust to fit the 
data (see Figure 9 below). 

The shapes of the fits in Figure 7 do not depend on the 
assumed inclination angle i. There is thus a two-parameter 
family of plausible models for the luminosity density, with 
a and i as free parameters. 

5.3 Determining the dust distribution 

Given a model for the luminosity density j, the optical depth 
of the dust can be determined with the following algorithm. 

Let the ra-axis be along the line-of-sight and let the z- 
axis (the symmetry axis of the galaxy) be oriented such that 
the observer sees the region z > in front of the dust disc 
and the region z < behind the dust disc. Define 



J 



jdw, I-{x,y) 



J 



j dw, 



(5) 



where (a;,?/) are Cartesian coordinates on the sky. The pro- 
jected intensity is 

I{x,y) = I+{x,y) + I.{x,y)e-^^'''y\ (6) 

where r(a;,2/) is the optical depth of the dust as seen pro- 
jected onto the sky. The observed intensity is obtained by 
convolving with the normalized seeing PSF, Plr): 



IobB{x,y) = J I{x' ,y')P{\f' - r\) dx' dy' , 



(7) 



where r = {x,y) and r" = [x',y'). 

For any given inclination and model for j one can eval- 
uate /-I- and /_, whereas Jobs and P(7') are known from 
observations (Section 2). To recover the unknown r(a;,2/) 
we define and evaluate 

I{x,y) = Iobs{x,y) - j I^{x' ,y')P{\T' - t\) dx' dy' . (8) 
Combination of equations (6)-(8) yields 

I{x,y) = J T{x,y)P{\r' -r\)dx' dy', (9) 
where we defined 

T(^,2/) = e-^(='^'/_(^,2/). (10) 

Recovering the unknown T(a;,2/), and hence r(a;,2/), from 
equation (9) is a seeing deconvolution problem. To solve 
it, we implemented a version of the Richardson-Lucy algo- 
rithm (Lucy 1974). Convergence was usually achieved after 
~ 40 iterations. Results are discussed in Section 5.5 (see 
Figure 8). 

5.4 Modelling the ionized gas distribution 

The observed ionized gas distribution on the plane of the sky 
determines the inclination angle i of our models, because the 
intrinsic distribution is assumed to be circularly symmetric 
in the equatorial plane. Recovering i from Figure 3 is not 
straightforward though, because: (i) the observed distribu- 
tion has been convolved with a circular seeing PSF; (ii) part 
of the line-emission flux might have been absorbed by dust; 
(iii) the transmission of the emission line flux by the on-band 
filter used in the narrow-band imaging varies with position, 
due to the rotation of the gas; and (iv) the flux distributions 
for Ha and the [Nil] lines are different. To circumvent prob- 
lems (iii) and (iv), we determine i by fitting to the fluxes ob- 
served in our major and minor axis spectra (Figure 6), rather 
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Figure 9. The dust optical deptli at tlie projected position of the 
galaxy centre, Tcenj as function of the cusp-steepness ck. Models 
with a steeper power-law cusp (smaller ck) imply a larger optical 
depth for the dust. 

than to our narrow-band image (Figure 3). The following 
approach was adopted to take (i) and (ii) into account. 

Consider a fixed choice for a, and a trial value for i. The 
parameters jo, b and /3 of the luminosity density are fitted 
as in Section 5.2. Subsequently, the optical depth r(a;,2/) 
is determined as in Section 5.3. A function f (-R) of the 
form (4) is then fitted to the integrated line flux observed 
in the major axis spectrum, taking into account the seeing 
convolution, the slit width and the binning into pixels along 
the slit. Given the fitted form for f (-R), the predicted fluxes 
are calculated for the minor axis spectrum, and compared 
to the data. This procedure is repeated for different values 
of i. The one that yields the best fit to the minor axis flux 
distribution is adopted as final estimate. 

5.5 Results 

We calculated the best-fitting inclination for Ha and [Nil] 
separately, and for both assumptions for the absorption of 
the emission-line flux by dust, fobs = F and fobs = F 
(as discussed in Section 5.1). In addition, we examined dif- 
ferent values of a (which comes in because it determines r). 
The best fit (see Figure 10 below) was without exception 
i = 70° ± 5°. 

Having fixed the inclination, T{x,y) depends on a only. 
Figure 8 shows the result of calculating it with the algo- 
rithm of Section 5.3, for two representative cases: a = and 
a = —1.3. For a = the optical depth is very asymmet- 
rical, the point of maximum r being offset from the centre 
of the outer isophotes by ~0.7" along the major axis, and 
~0.3" along the minor axis. The optical depth distribution 
for a = —1.3 is more symmetric (though not entirely), and 
the r-contours are more nearly elliptical. Assuming the dust 
distribution in the equatorial plane to be circularly symmet- 
ric, the outermost r isophote indicates i Ki 72°, close to the 
value derived from the ionized gas distribution. This was 
not imposed a priori, and hence provides a successful con- 
sistency check on the model. 

Figure 9 shows the required optical depth in the centre. 
Teen, as fuuctiou of the assumed cusp steepness a. Models 
with smaller a have more light near the centre, and thus 



require a larger optical depth for the dust to fit the data. 
Unfortunately, we only have broad-band images in one band. 
Independent constraints on the optical depth of the dust 
(and hence on a) from the observed reddening are thus not 
available. The observed infrared emission from NGC 7052 
(part of which might in fact be non-thermal) can be used to 
estimate the total mass of dust in NGC 7052 (Jura 1986), 
but this does not yield useful constraints on a. 

HST observations have revealed that essentially all el- 
liptical galaxies have central power-law cusps, rather than 
constant surface brightness cores (Crane et al. 1993; Fer- 
rarese et al. 1994). This favours models with a < —1 as 
most plausible explanation of our data. Such models also 
imply a more symmetrical dust distribution than models 
with, e.g., a = 0. In the following we restrict ourselves to 
the model with a = —1.3 (which, in the absence of dust 
absorption, corresponds to / oc r~°'^ in projection). As will 
be shown (see Figure 12 below), this yields the best fit to 
the observed major axis gas rotation curve. In addition, it 
is a more-or-less typical value for bright elliptical galaxies 
(Kormendy et al. 1994). For example, M87 has a = -1.26 
(Lauer et al. 1992). 

With i = 70° and a = —1.3 the best fitting values for 
jo, b and /3 are: jo = 6.0Lqpc~^, b = 2.35" and /3 = —0.58. 
We henceforth refer to this model as our 'standard model'. 
Its total luminosity is formally infinite. This is no great 
drawback. The agreement between the model and the data 
only breaks down at large radii, where no kinematic data 
are available anyway. 

Figure 10 shows the best fits to the integrated Ha fluxes 
along the major and minor axes for our standard model, with 
r(a;,2/) as in Figure 8, and f (r) as in equation (4). Models 
that assume fobs = f and fobs = e~^f are both shown. 
The latter assumption provides the better fit. It reasonably 
well explains the slight asymmetries in the flux distribution, 
as due to absorption by the somewhat asymmetric dust dis- 
tribution. In the following we restrict ourselves to this case. 
However, adopting the alternative assumption, fobs = f , 
would not change the main conclusions of our paper. 

The (fobs = e~^f ) model for the Ha fluxes in Figure 10 
has Ri = 0.1", R2 = 1.0", R = 1.5" and h/h = 120. The 
ratio of the total flux in the first exponential component of 
F{R), fi, to that in the second component, F2, is F1/F2 = 
2.7. Hence one can think of the flux as being the sum of 
a spatially unresolved component [Ri <C seeing FWHM) 
which contains most of the flux, and a resolved extended 
component. The observed flux in the [Nil] lines is more 
centrally peaked than that in Ha (Figure 6). It can be fit 
with a similar flux distribution as above, but with F1/F2 = 
13.5. The spatially unresolved component might be a so- 
called 'narrow-line region'. Such regions are seen in many 
active galactic nuclei, and have a typical size of ~100pc 
(e.g., Osterbrock 1991). This would indeed be unresolved at 
the distance of NGC 7052 (100 pc = 0.44"). 

Because the observed flux distribution can be reason- 
ably well fit by a model in which the intrinsic flux distri- 
bution is circularly symmetric in the equatorial plane, we 
have not studied any models in which the intrinsic distri- 
bution is not circularly symmetric. Such models are not 
implausible though, and could probably fit the data equally 
well. If photoionization is the ionization mechanism in the 
nucleus of NGC 7052, shielding of the ionization source by 
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Figure 10. The Hck emission- line flux observed in our major and minor axis spectra, as function of galactocentric distance. The vertical 
scale is in arbitrary units. The curves are the best (seeing convolved and spatially binned) flts of a circularly symmetric flux distribution 
in the equatorial plane, with i^(-R) as in equation (4) and i = 70°. The dashed curve assumes there is no absorption of line emission 
by dust. The solid curve does take into account the absorption of line emission by dust (i^obs = > with r as in the right panel of 

Figure 8). The latter model reasonably well explains the asymmetries in the observed flux distribution as being due to dust absorption. 



the slightly asymmetric dust distribution would lead to a 
slightly asymmetric flux distribution. The flux would be 
highest on the side with the smallest dust optical depth, in 
agreement with the observations. 



6 DYNAMICAL MODELS 

6.1 Modelling the stellar kinematics 

The stellar velocity dispersion determined in Section 4.1 con- 
strains the mass-to-light ratio of the stellar population. 

We assume the stellar phase space distribution function 
to be of the particular form f{E, Lz), where E is the energy 
and Lz the angular momentum around the symmetry axis. 
The predicted stellar kinematics can then be calculated from 
the Jeans equations. We did this using the technique and 
software described in van der Marel et al. (1994b), taking 
the seeing, slit width and spatial binning of the data into 
account. Our standard model requires a mass-to-light ratio 
[M/L^R = 5.0 ± 1.0 (in solar units) to reproduce the ob- 
served velocity dispersion of 266 ± 26kms~'" between —5" 
and 5". This value falls on the relation between mass-to- 
light ratio and luminosity obtained by van der Marel (1991) 
for a sample of 37 bright elliptical galaxies. The inferred 
mass-to-light ratio depends only little on either a or the 
possible presence of a central black hole, since these influ- 
ence only the dynamics of the stars very close to the nucleus 
(e.g., see Figure 11 below). 

Our data do not constrain the stellar dynamical struc- 
ture of NGC 7052, because we have not measured the vari- 
ation of the stellar kinematics as function of position on the 
sky. Our mass-to-light ratio estimate can be slightly in error 
if the dynamical structure differs from that of an f{E, Lz) 
model. In general, three-integral or triaxial models are more 
appropriate. However, f[E,Lz) models are much easier to 
construct, often provide reasonable (though not perfect) ap- 
proximations to the dynamics of real galaxies (van der Marel 
1991), and at least take proper account of flattening. 



6.2 Modelling the gas kinematics 

The circular velocity in the equatorial plane, Vc{R), satisfies 

K=(i^) = i^g, (11) 

where $ is the gravitational potential of the system. We 
assume that the potential consists of contributions from the 
stellar mass density ($*) and a possible central black hole: 
$ = $t — GM-bb.It. The total mass of the ionized gas in 
elliptical galaxies is typically of the order of 10^ - 10^ Mq 
(Phillips et al. 1986), so we neglect its contribution to the po- 
tential. The stellar mass density in our models, p = [M/L^j, 
is spheroidal with constant flattening. The components of 
the force V$t exerted by such a mass density can be ex- 
pressed as one-dimensional quadratures (Binney & Tremaine 
1987, equation [2-88]) that are easily evaluated numerically. 
There is thus no need to make the simplifying assumption 
that the potential is spherical, as is often done (e.g., Cald- 
well, Kirshner & Richstone 1986). 

In an axisymmetric potential cold gas will eventually 
settle in the equatorial plane (e.g., Tohline, Simonson & 
Caldwell 1982). In triaxial potentials the case is more com- 
plicated (e.g., de Zeeuw & Franx 1989). Near the nucleus 
the settling time is generally much shorter than the Hubble 
time, although this need not be the case in the outer parts 
(see de Zeeuw 1994 for a review). Here we assume that NGC 
7052 is axisymmetric and the gas has indeed settled and 
reached equilibrium. Hence, it will move on circular orbits 
with velocity Vc{R). Figure 11 shows the line-of-sight ve- 
locity component ■uios = V^(i?)sini, as function of projected 
distance along the major axis, for models with several differ- 
ent black hole masses. If present, the black hole dominates 
the potential close to the centre and induces a Keplerian rise 
in the circular velocity, V^^R) oc R~^^^ for R ^ 0. 

In addition to its gravitational motion, the gas has a 
(one-dimensional) thermal velocity dispersion 
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Figure 11. The mean line-of-sight velocity of tlie gas along 
the major axis, uios = Vc{R)sin.i (in kms~^), for our standard 
model. The [M/L^ of the stellar population is fixed as deter- 
mined in Section 6.1. Predictions are shown for several different 
masses of a possible nuclear black hole : Afgn = 0, 1.5 X 10® , 
5 X 10* M0 and 1.5 X 10^ Mq. The abscissa is in pc, with l" on 
the sky indicated by an arrow. 

where k is Boltzmann's constant, and m is the average mass 
of a gas particle. The typical temperature of the ionized 
gas-component in LINERS is ~10^ K (Osterbrock 1989), 
and iTthermai is thus On the order of ~10kms~''. 

The gas might also have a non-thermal velocity disper- 
sion. We will loosely refer to such a contribution as 'tur- 
bulent', although we will not attempt a physical description 
of the motions. We merely assume any turbulence to be 
isotropic, and parametrize the (one-dimensional) turbulent 
velocity dispersion phenomenologically as 



fturb(-R) = iJo X exp[ — R/Rt), 



(13) 



where Rt is the scale length of the turbulence. Turbulent 
motions could arise through, e.g., shocks or interaction of 
the gas with the radio jet of NGC 7052. Note in this context 
that, although LINERS are generally thought to be ionized 
by a non-thermal continuum, shocks have often been dis- 
cussed as alternative or additional source of ionization. See 
Ho, Filippenko & Sargent (1993) and Maoz et al. (1995) for 
recent reviews and discussions on this subject. 

The thermal and turbulent velocity dispersion need 
not disturb the bulk flow of the gas on circular orbits. 
We therefore assume the local velocity distribution of the 
gas, C{R, v), to be Gaussian with mean line-of-sight veloc- 
ity ■uios = V^(i?)sini, and line-of-sight velocity dispersion 



C{R,v) 



exp 



K(J?)i 



(14) 



The gas could also have a non-zero gravitational (bulk) 
dispersion. This would be the case if, e.g., the gas resided 
in individual clouds orbiting on non-circular orbits in the 
gravitational potential of the system (as suggested by, e.g., 
Cinzano & van der Marel 1994, in their models for NGC 
2974). Although physically quite different from turbulence. 



it leads to similar predictions for the observed gas kinemat- 
ics. The only difference is that the predicted mean streaming 
velocities are somewhat lower, due to the effect of 'asymmet- 
ric drift' (e.g., Binney & Tremaine 1987). We show below 
that models with turbulent velocity dispersions fit the ob- 
served rotation curve quite well. In the following, models 
with non-zero bulk velocity dispersion for the gas are there- 
fore not discussed. 

For any particular position along the slit we adopt a 
Monte-Carlo procedure to take into account the effects of 
seeing, the slit width and pixel size. At each Monte-Carlo 
step a random position is drawn from the rectangular re- 
gion on the sky defined by the slit width and pixel size, and 
a two-dimensional vector is added, drawn at random from 
the PSF. The line-of-sight velocity distribution at this po- 
sition (eq. [14]) is then calculated, and added to the sum 
of the previous Monte-Carlo steps, weighted with the flux 
(taking dust absorption into account). Some 10^ steps suf- 
fice to obtain accurate results. The result is convolved with 
a Gaussian of the appropriate width to model the instru- 
mental velocity resolution. This then yields the predicted 
emission line profile shape as function of velocity. 

As a brief aside, let us stress here the importance of 
modelling the seeing. Consider for example the case of a 
galaxy with a homogeneous mass density core, and with 
constant M/L. The slope of the circular velocity curve is 
then proportional to the core density. Seeing flattens the 
slope of the observed rotation curve near the nucleus. So 
if the core density is determined from the observed slope 
without modelling the seeing, it is underestimated. As a 
result, it will be inferred incorrectly that M/L has a central 
minimum. This effect might be responsible for the increase 
in M/L with radius reported by Caldwell et al. (1986) for 
NGC 7097. This result was previously questioned by Bertola 
et al. (1991), who suggested it to be due to gas motions on 
elliptic orbits in a triaxial potential. 

6.3 Gas rotation velocities and velocity dispersions 

Figure 12 shows the rotation velocities and velocity disper- 
sions predicted by our model, for Ha and [Nil] separately, 
along both the major and minor axes. Solid curves are for 
the case when there is no central black hole and no gas tur- 
bulence. The major axis rotation curves are extremely well 
fit. The predicted slope of the rotation curve depends on the 
stellar mass near the nucleus, which is determined mainly 
by a. The value a = —1.3 yields the best fit, which is why 
we adopted this value for our standard model. The fit is less 
good for other values of a, although one could improve the 
fits of such models by choosing a slightly different value for 
the stellar mass-to-light ratio (assuming that the estimate 
in Section 6.1 would not be correct). Either way, no black 
hole is required to fit the steep major axis rotation curve. 
However, the model fails to explain the central peak in the 
observed gas velocity dispersions. This was found to be the 
case independent of the particular choice for a or the mass- 
to-light ratio. To fit the observed velocity dispersions of the 
gas one needs to invoke the presence of a massive black hole, 
gas turbulence, or a combination of both. 

Seeing convolution of a Keplerian rotation curve around 
a massive black hole increases the observed central velocity 
dispersion. Dashed curves in Figure 12 show the velocity 
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Figure 12. Rotation velocities and velocity dispersions for Hck and [Nil], along the major and minor axes. The solid curves are the 
predictions of our standard model, without gas turbulence and without a central black hole. The major axis rotation curves are well fit. 
The predicted minor axis rotation velocities are slightly non-zero because of the dust absorption in the model. The model does not fit 
the observed central peak in the major axis velocity dispersion. This can be remedied by invoking gas turbulence or the presence of a 
central black hole. The dashed curves in the bottom panels show the effect of adding a black hole of mass 1.5 X 10^ , 5 X 10^ and 
1.5 X 10^ M0 , respectively. The central velocity dispersion for [Nil] is higher than for Hck, because the flux distribution of the former is 
more centrally concentrated. Inclusion of a black hole in the model does not alter the predicted rotation velocities signiflcantly. 



dispersions for NGC 7052 predicted by models with different 
black hole masses. The observed dispersions are reasonably 
well fit when Mbh = 10^ Mq. Inclusion of a black hole 
does not alter the predicted rotation velocities significantly, 
because the Keplerian part of the circular velocity curve is 
much smaller than the seeing radius (cf. Figure 11). 

Alternatively, one can fit the observed velocity disper- 
sions by invoking gas turbulence of the form (13). The ob- 
served width of the central peak in the gas velocity disper- 
sions sets an upper limit on the scale length of the turbulence 
Rt- Good fits to the data can be obtained for all values of 
Rt smaller than the seeing FWHM of the major axis ob- 
servations. If Rt = 0.5", one needs a central turbulence 
iTo ~ 350 km s"''. If Rt is equal to the scale length of the 
unresolved part of the flux distribution (0.1", see Section 
5.5), one needs ao ~ 1000 km s"''. 

To examine whether turbulent motions on the order of 
1000 km s"'" are realistic, we calculated the central escape 
velocity for our standard model without a black hole. A gen- 
eral expression for this quantity for models with luminosity 
densities of the form (2) is given by Qian et al. (1995). The 
central escape velocity is 1459 km s"'^. It is finite in spite of 
the infinite mass of our model. Turbulent motions exceeding 
or even close to the central escape velocity are unlikely to 
exist. This argument favours models with smaller ao and 
larger Rt. 

The models with either a central black hole or turbu- 
lence both predict a higher central velocity dispersion for 



[Nil] than for Ha, as is observed. This is because the flux 
distribution of [Nil] is more peaked towards the nucleus than 
that of Ha. Consequently, the observed [Nil] flux in the cen- 
tral pixel originates on average from a region nearer to the 
nucleus (where the RMS velocity is higher) than the ob- 
served Ha flux. 

6.4 Emission line profiles 

The observed gas velocity dispersions can be fit either by 
invoking the presence of a central black hole, or by invok- 
ing gas turbulence. In an attempt to discriminate between 
these models we now consider the predicted and the ob- 
served emission line profile shapes. The observed line pro- 
files away from the centre are slightly asymmetric with a 
tail at low rotation velocities (Figure 5). This is caused by 
seeing convolution, and is reasonably well reproduced by all 
models that fit the major axis rotation curve. Of more in- 
terest here is the shape of the emission line profile in the 
very centre. 

Figure 13 compares model predictions (thick curves) 
and observations (thin curves) for the observed central Ha 
emission line. The left panel shows the predictions for the 
model with a 10^ Mq black hole. The predicted emission line 
shape is completely different from the observed line shape 
and we hence rule out this model. The predicted line shape 
has three local maxima. These can be understood in terms 
of the circular velocity curve of the model (Figure 11). Due 
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Figure 13. The thin curve in each panel is the observed Hck emission Une (as in Figure 5) for the galaxy centre (r = —0.05"). The thick 
curves show the predictions of our models. Left: Afgn = 10^ ^0) and no gas turbulence; Middle: gas turbulence with (Tq = 350kms~^ 
and Rt = 0.5", and no central black hole; Right: Gas turbulence with (Tq = 300kms~^ and Rt = 0.5", and also a 5 X 10^ central 
black hole. The vertical scale is in arbitrary units. The latter two models both provide an acceptable fit. 



to the Keplerian rise at small radii, there is no gas along the 
major axis with I'UiosI < 'Wmin ~ 300 km s"''. This results in 
two local maxima in the predicted line shape at velocities 
~±'Umin. The predicted flux at these two local maxima is not 
the same in this particular case, because: (i) the observed 
profile is slightly offset from the nucleus [r = —0.05"); and 
(ii) the model takes into account absorption by dust, which 
is asymmetrically distributed. The third local maximum in 
the line profile, at i; Ri Okms"'', is due to the gas outside 
the region influenced by the black hole, which is roughly in 
solid body rotation. If the black hole mass were reduced to 
zero, this would be the only remaining maximum. 

The middle panel of Figure 13 shows the predicted line 
shape for the model with a turbulent gas dispersion with 
iTo = 350 km s"'^ and Rt = 0.5". The shape is smooth and 
close to Gaussian, and fits the observed shape well. Hence, 
this model is consistent with the data. 

The left panel of Figure 13 indicates that any viable 
model with a black hole must have gas at velocities smaller 
than ■Umin (as defined above), to avoid the prediction of local 
maxima in the line profile shape (which are not observed). 
This can be achieved by constructing combined models, in 
which the gas near the nucleus is in Keplerian motion around 
a black hole, but in addition has turbulence of order ■Umin. 
We have constructed several of these combined models. An 
adequate fit to the observed velocity dispersions can be ob- 
tained with, e.g., Mbh = 5 x 10* Mq and turbulence with 
ITo = 300 km s"'^ and Rt = 0.5". The right panel of Figure 13 
shows the central line profile predicted by this model. It fits 
the data as well as the model in the middle panel, which has 
turbulence but no black hole. 

6.5 Predictions for HST observations 

Our ground-based observations show that the motions of 
the ionized gas in NGC 7052 cannot be purely circular, but 
must have a significant random (turbulent) component near 
the centre. A central black hole of mass Mbh ^ 5 x 10* Mq 
might or might not be present. 

High spatial resolution spectra with the HST could pro- 



vide more stringent constraints on the presence of a black 
hole. We calculated predicted rotation velocities and veloc- 
ity dispersions for hypothetical observations with the 0.26" 
diameter circular aperture of the Faint Object Spectrograph 
(FOS) on the HST. Figure 14 shows the results for the mod- 
els discussed in Sections 6.3 and 6.4, as function of posi- 
tion along the major axis. Models with a black hole predict 
higher rotation velocities than those without a black hole. 
One does not expect to observe rotation velocities as high as 
have been observed in M87 (Harms et al. 1994), because M87 
is closer by than NGC 7052, and has Mbh w 2.4 x 10^ Mq. 

7 CONCLUSIONS 

We have presented high spatial resolution broad-band imag- 
ing, narrow-band imaging and spectroscopy for the central 
region of the E4 radio galaxy NGC 7052. The narrow-band 
image and spectra show that the (previously unresolved) 
ionized gas near the nucleus resides in a rapidly rotating 
disc with radius ~1.5". The emission line spectrum is char- 
acteristic of that of a LINER. We have determined the kine- 
matic quantities and emission line shapes of the gas along 
both the major and minor axes. The observed major axis 
rotation velocity rises to ~250kms~'" at ?• = 1.5". The ve- 
locity dispersion of the ionized gas is ~210 — 260 km s"'' in 
the centre and ~60kms~'' at ?• = 1.5". The central veloc- 
ity dispersion for [Nil] is higher than for Ha, because the 
flux distribution of [Nil] is more centrally concentrated than 
that of Ha. Our broad-band image confirms the previously 
reported presence of a nuclear dust lane, extending from 
r Ri —1.6" to 1.6" along the major axis. 

We have constructed detailed models to interpret the 
data. The galaxy is assumed to be axisymmetric with con- 
stant flattening, with the gas and dust in the equatorial 
plane. The gas distribution is assumed to be circularly sym- 
metric. We allow for the possibility of absorption of the 
emission-line flux by the dust. The optical depth distribu- 
tion of the dust on the sky, r(a;,2/), cannot be inferred di- 
rectly from the observed broad-band image, since it depends 
on the (unknown) fractions of the observed intensity that 
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Figure 14. Predicted rotation velocities and velocity dispersions along the major axis for observations with the 0.26" diameter circular 
aperture of the FOS on the HST. The three models of Figure 13 are shown, as well as the model with no black hole and no turbulence. 
Models with a black hole discern themselves by predicting larger rotation velocities. HST observations can thus provide more stringent 
constraints on the presence of a black hole in NGC 7052. 



originated in front of and behind the dust 'screen', respec- 
tively. However, t{x, y) is fixed uniquely by the observations 
once a model for the three-dimensional light distribution has 
been specified. We present a new algorithm for recovering 
r(a;,2/), that deconvolves simultaneously for the effects of 
seeing. 

The observations of the ionized gas disc, when corrected 
for the effects of seeing, indicate an inclination angle i Ki 70°. 
The flux distribution of the ionized gas is extremely centrally 
concentrated. The flux inside the seeing radius is unresolved, 
hinting at the presence of a narrow-line region in this galaxy. 
The central cusp-steepness a of the stellar light distribution 
[j oc for ?• ^ 0) is not well constrained by our data, 
because of the absorption by dust near the centre. Moti- 
vated by a variety of arguments we choose a = —1.3 in our 
'standard model'. The observed dust morphology is then 
consistent with a roughly circularly symmetric dust distri- 
bution in the equatorial plane, with maximum optical depth 
in the centre; i.e., a disc, rather than a ring or filament. 

To model the gas kinematics, we assume the gas to move 
on circular orbits in the equatorial plane. The local velocity 
dispersion of the gas is assumed to be the sum of a thermal 
and a turbulent (or otherwise non-gravitational) component. 
The circular velocity is calculated from the combined grav- 
itational potential of the stars and a possible nuclear black 
hole. The average mass-to-light ratio of the stellar popula- 
tion is assumed to be constant. Its value is determined by 
fitting to the stellar velocity dispersion determined from our 
data, assuming a stellar phase-space distribution function of 
the form f{E, L^). 

The model with no black hole and no turbulence fits the 
major axis gas rotation curve, but fails to fit the observed 
central peak in the gas velocity dispersion. This indicates 
that there must be a strong increase in the RMS gas motion 
inside the seeing radius. A model with a 10^ Mq black hole 
and no turbulence yields a good fit to the observed veloc- 



ity dispersions, due to seeing convolution of the Keplerian 
part of the rotation curve. However, this model is ruled out 
because it predicts a nuclear emission line profile that has 
two pronounced peaks at i; Ri ±300 km s"'', which are not 
observed. Models without a black hole, but with strong tur- 
bulence (> 300 km s"'') within the seeing radius, can fit not 
only the rotation curve and velocity dispersions, but also the 
observed line profile shapes. Combined models with both a 
black hole and strong turbulence can also fit the data, but 
the black hole mass must be < 5 x 10* Mq. Predictions for 
future observations with the HST show that such observa- 
tions should yield more stringent constraints on the mass of 
a possible black hole. Any black hole with mass > 10* Mq 
will reveal itself by causing the rotation curve observed with 
HST to be significantly steeper than that reported here. 

A complication not addressed here is the possibility of 
elliptic gas orbits in a triaxial potential. Ideally one would 
like to include this in the modelling. This has in fact been 
done successfully for extended ionized gas disks in several 
early-type galaxies (Bertola et al. 1991; Amico et al. 1993). 
However, information on the two-dimensional velocity field 
on the sky is required to place meaningful constraints on 
such models. Such information is not available for the galaxy 
we study here, and will generally be very difficult to obtain 
for any gas disk restricted only to the nuclear region. 

Although the upper limit we have set here on the mass 
of a black hole in NGC 7052 is not very stringent, it is al- 
ready ~5 times smaller than the black hole mass inferred for 
M87 from HST data (Harms et al. 1994). It will be worth- 
while to search for and observe other galaxies with nuclear 
gas discs that are closer by than NGC 7052. Such galaxies 
can be studied at higher spatial resolution, as measured in 
physical units, allowing the detection of less massive black 
holes. One galaxy for which HST observations are planned 
is NGC 4261. HST photometry of this galaxy in the Virgo 
cluster have revealed a beautiful nuclear dust disc (Jaffe et 
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al. 1993). Emission line spectra obtained with the WHT 
show a broad emission-line component, which hints at the 
presence of a ~4 x lO'^ Mq black hole. The observed nu- 
clear emission line shapes are not double (or triple) peaked, 
so some form of gas turbulence is again required (Jaffe et 
al. 1994). In the near future, kinematic observations of nu- 
clear gas discs are likely to become a widely used tool in the 
search for massive black holes in galactic nuclei. The mod- 
elling and analysis techniques developed here will be useful 
for interpreting such data. 
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